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ABSTRACT 

We have determined the spectral energy distribution at wavelengths between 6 cm and 
850 fim for the proto- typical S(stellar)-type symbiotic star, CI Cygni, during quies- 
cence. Data were obtained simultaneously with the Very Large Array and the SCUBA 
sub-millimetre (sub-mm) camera on the James Clerk Maxwell Telescope. The data 
have allowed us to determine the free-free turnover frequency of the ionised compo- 
nent, facilitating a model-dependent estimate of the binary separation to compare 
with the known orbital parameters of CI Cyg and to critically test the known models 
for radio emission from symbiotic stars. In particular, our data rule out the two most 
popular models: ionization of the giant wind by Lyman continuum photons from its 
hot companion and emission resulting from the interaction of winds from the two bi- 
nary components. 
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1 INTRODUCTION 

The symbiotic stars are a collection of objects which exhibit 
the characteristic optical-infrared absorption bands of a cool 
giant upon which are superimposed emission lines whose 
origin requires conditions typical of a much hotter object. 
Multi-frequency observations have led to the universal ac- 
ceptance that they are binary systems, typically containing 
a post-asymptotic giant branch star or main-sequence dwarf 
accreting matter from a late- type giant or supergiant. 

Radio emission from symbiotic stars is dominated by 
free-free radiation from gas ionised by the hotter of the bi- 
nary companions. Seaquist, Taylor & Button (1984, here- 
after STB) and Seaquist et al. (1990, 1993) investigated 
the radio properties of ~ 100 symbiotics, finding differences 
between the ionised regions around D-type systems (those 
containing Miras, where the IR emission is predominantly 
from dust) and S-type systems (containing first-ascent gi- 
ants where the IR emission is predominantly photospheric) . 
The SEDs of D-type systems turn over from partially opti- 
cally thick (a ~ 1, where F v oc v a ) to optically thin emission 
(a ~ —0.1) at centimetre wavelengths, whereas the SEDs of 
S-type systems remain optically thick throughout the radio 
regime. 

The frequency of the turnover to optically thin free-free 
emission, f t , is of particular interest: in the STB model, i/ t oc 
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a~ 7 (where a is the binary separation) and the optically 
thin portion of the free- free spectrum is very flat (a ~ 0); 
in the colliding winds model, u t oc a -1,5 , and the spectral 
index never flattens (a ~ +0.5, even for v > v t ). 

Mm/sub-mm observations have determined the spec- 
tral turnovers in several D-type systems (Ivison, Hughes & 
Bode 1992; Seaquist & Taylor 1992; Ivison et al. 1995). Until 
recently, however, the limited sensitivity of bolometer de- 
tectors had hindered attempts to pin down the turnover 
frequency in quiescent S-type systems. These are thought 
to be fainter, probably because of a feedback between the 
lower mass-loss rate of their earlier-type giant (compared 
with D-type giants) and the consequent low accretion rate 
onto the ionising companion star. The result is a lower gas 
density and a lower flux of ionising radiation, which both 
lead to a lower emission measure. Brighter S-type systems 
(e.g. AG Peg, PU Vul) are often recovering from nova-type 
outbursts making them unsuitable representives of the nor- 
mally quiescent population. 

This situation is unfortunate because the STB model 
can be tested against the SEDs of quiescent systems: a crit- 
ical parameter in the model is the binary separation, which 
is known for many S-type symbiotic binaries (where periods 
are a few years, typically) but not known for D-type systems 
(which are thought to be more widely separated and to have 
periods spanning one or more decades). 

Attempts have been made to constrain i/ t for several 
well-studied S-type systems (RW Hya — Mikolajewska & 
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Table 1. Observed properties of CI Cyg. 



Property 

a(J2000) 
<5(J2000) 

3600 A (17) 
4400 A (B) 
5500 A (V) 
7000 A (R) 
1.25 urn (J) 
1.65 /im (H) 
2.20 /um (A") 

12 /im 
25 /im 
60 fim 
100 ^trn 

450 fim 
850 

1350 (im 
2000 (im 

0.69 cm 
1.33 cm 
2.01 cm 
3.54 cm 
6.17 cm 



Telescope 
VLA 

CAO 



UT Date 
1998 Feb 15 

1985 Jun 07 



IRAS 



JCMT 



VLA 



1984 



1998 Jan 29, 31 



1998 Jan 30 



1998 Feb 20 



1998 Feb 15 



Measurement 

19 h 50 m 11.83 s 
+35°41'03.2" 

16.8 ± 0.1 mjy 
44.4 ± 0.1 mjy 
0.16 ± 0.01 Jy 
0.77 ± 0.01 Jy 
8.02 ±0.07 Jy 
10.4 ± 0.1 Jy 
10.1 ± 0.1 Jy 

0.93 ±0.06 Jy 
0.33 ±0.03 Jy 
3(7 < 0.45 Jy 
3cr < 5.97 Jy 

3cr < 42.7 mjy 
9.50 ± 1.29 mjy 
10.1 ± 1.4mJy 
11.0 ± 2.8 mjy 
7.35 ± 1.55 mjy 
3.92 ± 0.75 mjy 
3.01 ± 0.29 mjy 
1.49 ± 0.10 mjy 
0.91 ± 0.08 mjy 



Comment 

From 3.54-cm image, ±0.1". 

East of position reported by Seaquist et al. (1993). 
Munari et al. (1992). 



Whitelock & Munari (1992). 



Photometry mode. 



From raw 0-150 kA visibilities. 



From map made with 0-150 k A visibilities. 



Omont 1998; AG Peg and Z And — Ivison et al. 1992, 1995) 
but these studies have usually had to rely on estimates of 
the optically thin free- free emission in the sub- mm range de- 
rived from optical/ultraviolet (UV) Hi free-free and bound- 
free emission measures. The consistency between the data 
and binary models revealed by the aforementioned work is 
not therefore based on direct measurements of the turnover 
frequency, only on an optically thin emission measure in- 
ferred from optical/UV spectroscopy. We need to sample 
the mm/sub-mm spectral energy distribution of S-type sym- 
biotic stars in order to define the turnover frequency in a 
model-independent manner. 

For the investigation reported here, we selected the well- 
studied S-type system, CI Cyg, a proto-typical symbiotic 
star with a known orbital period (P or b = 855.6 d) and a 
well-defined spectroscopic orbit (Mikolajewska 1997; Bel- 
czyriski et al. 2000 and references therein). A thorough spec- 
troscopic study of CI Cyg by Kenyon et al. (1991) demon- 
strated that it consists of an M5II asymptotic branch giant, 
M g ~ 1.5 Mq, and a ~ 0.5 Mq companion separated by 
2.2 AU. Our objective was to define the shape of the spec- 
tral energy distribution from centimetre (cm) to sub-mm 
wavelengths, determining v% and thus quantitatively testing 
and discriminating between the existing models for symbi- 
otic binaries. In §2, we report our measurements of CI Cyg, 
which we obtained near-simultaneously at wavelengths be- 
tween 850 pm and 6 cm during 1998 January and February. 
Our interpretation of the resulting spectral energy distribu- 
tion is presented in §3. 



2 OBSERVATIONS 



2.1 Measurements at 0.7 — 6.2 cm with the VLA 

Observations were carried out during 1998 February 15 and 
20 (ut) with the NRAO Very Large Array (VLA), New 
Mexico, during a move between the D and A configura- 
tions. Around 30min was spent on source at 0.69, 1.33, 
2.01, 3.54 and 6.17 cm, with measurements of the bright 
gain/phase calibrator, 2015±371, every 5-10 min. The flux 
density scale was set using observations of 3C 286 at 2.01, 
3.54 and 6.17 cm, and using a variety of bright calibrators 
at 0.69 and 1.33 cm. We recorded a bandwidth of 100 MHz 
- two contiguous 50-MHz bands, with both right and left 
circular polarizations. 

Calibration of the synthesis data followed standard 
NRAO Cookbook procedures within MPS. At 6.17 and 
3.54 cm, the MX routine was used to make and clean maps, 
restricting the uv coverage to 0-150 kA to reject those an- 
tennas which had already been moved into A-configuration 
positions. At shorter wavelengths we utilised prior knowl- 
edge of the source position (from a 3.54-cm map made using 
all baselines) to determine the flux density directly from the 
visibility data, again restricting ourselves to short baselines 
typical of D configuration. The measured flux densities are 
listed in Table 1. 



2.2 Measurements at 450 — 2000 ^m with SCUBA 
on the JCMT 

Observations with SCUBA (Holland et al. 1999) on the 15- 
m James Clerk Maxwell Telescope (JCMT), Mauna Kea, 
Hawaii, were made using its photometry mode at 450, 850, 
1350 and 2000 /jm during the nights of 1998 January 29- 
31 (ut). The secondary mirror was chopped in azimuth by 
60" at 7.8 Hz and jiggled every 1 s in a simple 9-point pat- 
tern with 2" offsets. The telescope was nodded between the 
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Figure 1. The spectral energy distribution of CI Cyg between 
the radio and optical wavebands. New measurements are repre- 
sented by empty circles, old measurements by filled squares. The 
dashed and the dotted lines, respectively, represent power laws 
with spectral indices +0.96 and — 0.10 (see text). 



signal and reference beams every 9 s in a signal-reference- 
reference-signal pattern. The pointing was checked regularly 
and skydips were performed every hour to measure the at- 
mospheric opacity. 

Reducing the data, spikes were carefully rejected then 
measurements in the reference beam were subtracted from 
those in the signal beam. The data were corrected for at- 
mospheric opacity and calibrated against Uranus and Mars. 
The measured flux densities are listed in Table 1. 



3 RESULTS AND DISCUSSION 

The radio spectrum of CI Cyg is shown in Table 1 and 
Fig. 1 together with earlier radio data (Torbett & Campbell 
1989; Seaquist, Krogulec & Taylor 1993) and published opti- 
cal/IR measurements (Munari et al. 1992; Whitelock & Mu- 
nari 1992). The resulting spectral energy distribution covers 
the entire radio to optical range. The IR region is clearly 
dominated by photospheric emission (from an M5 III giant 
- solid line in Fig. 1) while the cm to sub-mm emission is 
dominated by partially optically thick bremsstrahlung (free- 
free) from ionised gas. The radio spectrum turns over at mm 
wavelengths indicating that it has become optically thin to 
free- free emission. 

To estimate the turnover frequency, v t , the (partially) 
optically thick cm-wave radio data were fitted with a power 
law spectrum, Si, oc v a . The optically thin sub-mm data 
were similarly fitted but the spectral index was arbitrary 
fixed at a = —0.1 (appropriate for optically thin emis- 
sion). For v t , we have adopted the point at which the 
low-frequency extrapolation of the optically thin free-free 
omission exceeds the observed optically thick emission by 
a factor e (r„ = 1). A least-squares fit to our data yields 
S* hick = 0.20 ± 0.02(^/GHz)°' 96(±0 ' 03) mJy and S* hin = 
17.5 ± 0.6(^/GHz) -01 mJy for the cm and sub-mm data, 
respectively, and we estimate v t = 27 ± 5 GHz. 

For free-free emission in the optically thick regime, the 
effective radio photosphere emitting at frequency v and tem- 
perature T is given approximately by the blackbody radius, 



In the case of CI Cyg, the radius at the optically thick- 
thin transition (i.e. v = v t ) is 52.6 (T/10 4 k)~°' s AU, adopt- 
ing a distance of 2 kpc (Kenyon et al. 1991), which should 
represent the inner radius of its optically thick shell. This 
radius is significantly larger than the size of the known 
binary orbit of CI Cyg, for any reasonable temperature. 
The emission measure of the shell, derived from the op- 
tically thin turnover frequency, v% ~ 27 GHz, is around 
6 x 10 14 (T o /10 4 K) 1 ' 35 cm~ 6 AU measured along the line-of- 
sight through its densest region. 

Finally, the optically thin flux density is consistent 
with a volume emission measure of nlV w 3 to 5 x 
10 58 (d/2kpc) 2 cm~ 6 for the respective range of T c ~ 10 4 
to 10 5 K found in CI Cyg (Mikolajewska 1985; Kenyon et 
al. 1991). This emission measure is an order of magnitude 
lower than the emission measure of 3 to 10 x 10 59 cm -6 de- 
rived from the optical/UV bound- free and free- free contin- 
uum and Balmer Hi emission lines, adopting a reasonable 
range for E(B — V) (~0.3-0.45), thus the free-free radio 
emission and the hydrogen optical/UV recombination emis- 
sion from CI Cyg do not seem to arise from the same volume 
of gas. 

The radio emission also provides a model-independent 
lower limit to the Lyman continuum photon luminosity, 



> 7x io«^ (J^Y 0A (JLJ) - 1 (JL 

" 7X1U mJy ll0 4 J I GHz J ^kpc 



(2) 



(Seaquist & Taylor 1987), which gives Ll y 5 to 2 x 
10 45 phots _1 for the appropriate range of T c (10 4 to 
10 5 ), much lower than the Ll y ~ 10 47 phots _1 indi- 
cated by studies of the UV continuum and Hell emission 
lines observed during quiescence (1985-94; Mikolajewska & 
Holowacz 2000). 

The current radio spectrum may be used to test the 
known models for radio emission from symbiotic stars. 
Above all, the very flat shape of the optically thin portion of 
our spectrum rules out the colliding wind model (in which 
the spectral index is never expected to flatten: a ~ +0.5, 
even for v > v% (Kenny 1995). 



3.1 Single-star models 

The very steep optically thick spectral index is also inconsis- 
tent with an isothermal, spherically symmetric steady flow 
(for which a = +0.6, Wright & Barlow 1975). There are 
several possible reasons for such a steep spectral index. For 
example, Seaquist & Taylor (1987) showed that in the case 
of a wind with a power-law temperature and density profile 
of the form T oc r~ n and n oc r~ p , respectively, the observed 
spectral index is given by 

4p - 6.2 - 0.6n 



2p- 



1.35n 



(3) 



In principle, then, a = +0.96 can be accounted for by a 
steady wind (p = 2) and a temperature gradient with a 
power- law index, n = 1.55, or by an isothermal wind with a 
density gradient, p = 2.51. 

Seaquist & Taylor (1987) noted that the optically thick 
portion of the radio spectrum would be steeper than that 
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expected for a steady, spherically symmetric flow if the ob- 
ject underwent episodic mass loss (i.e. if the flow was not 
steady) . The turnover of the spectrum at high frequencies is 
then due to the finite inner boundary of a shell-like emission 
region caused by the episodic mass loss. 

Although only a few attempts have been made to de- 
tect CI Cyg at radio wavelengths, it seems that the radio 
emission has been variable over at least last two decades (Ta- 
ble 2) and this variability could in principle be related to the 
large eruption that gave rise to an optical maximum in 1975. 
CI Cyg was not detected at 6 cm in 1982 February, whereas 
in 1987 October its 6-cm flux density was a factor of 4 — 5 
stronger than the 3<r upper limit from 1982 (Table 2) and a 
factor ~ 2 higher than our present measurement. It seems 
unlikely, however, that the radio-emitting gas was ejected in 
1975 because the delay between the maximum of the opti- 
cal/UV light curve and the onset of the radio emission is too 
long. For comparison, this delay was of order of days (UV) 
to months (optical) for the outburst of Z And (Fernandez- 
Castro et al. 1995). Also, the present characteristic radius at 
the optically thin-thick transition cannot represent the inner 
radius of an optically thick shell ejected over 20 years ago 
unless the expansion velocity of this shell is extremely low 
(^ 13kms _1 ). Such a low expansion velocity is, of course, 
compatible with a wind from the cool giant and a radio spec- 
trum of the form observed for CI Cyg can be produced by a 
discontinuity in that wind. However, we consider such a sce- 
nario to be unlikely. In particular, if the 1975 outburst were 
due to an episodic increase in the cool giant wind (as in the 
model proposed by Nussbaumer & Vogel 1987) changes in 
the optical and near-IR magnitudes of the giant should have 
been observed, contrary to the observations (Mikolajewska 
& Kenyon 1992). 

The respective absence and presence of radio emission 
in 1982 and 1987 may be related to the spectral changes 
in the optical/UV continuum and emission-line light curves 
observed in the 1980s (Mikolajewska 1985; Kenyon et al. 
1991). In particular, in 1975-83 (during the outburst and 
its decline) eclipses in UBV continuum and optical Hi, He I 
and He II emission lines were narrow, with a total duration of 
^ 0.16 P or b with well-defined eclipse contacts. Mikolajewska 
(1985) noted that the radii of the H + and He +2 regions in 
the orbital plane (as derived from the eclipse profiles) were a 
factor of 2 — 3 smaller than the respective spherical radii de- 
rived from the observed emission measures, indicating that 
the emission may arise in a non-spherical (perhaps bipolar) 
nebula or flow. Since 1984, the eclipses in both the con- 
tinuum and the emission lines (except for Hen, where the 
eclipse behaviour has remained practically unchanged) have 
become very broad and the light curves show continuous si- 
nusoidal variations suggesting that the bulk of the Balmer 
Hi and Hel emission is coming from a region between the 
binary components, possibly with a density concentration 
towards the cool giant and in the orbital plane. 

The transition from narrow eclipses to sinusoidal varia- 
tions was accompanied by an increase by a factor of ~ 2 in 
Balmer H I line and continuum emission, a decline in [O ill] 
and [Nem] emission and the appearance of very strong, 
highly ionised [Fevn], [Nev] and [Mgv] emission lines in 
the optical and UV bands. At the same time, both the 
eclipse behaviour and the intensity of Hell lines remained 
unchanged. The physical conditions in the nebular region 



Table 2. History of radio observations of CI Cyg 



Date 


Freqency 


Flux 


Reference 




(GHz) 


(mjy) 




1982 Feb 9 


4.9 


< 0.42 


ST90 


1987 Oct 


4.9 


1.92 ±0.04 


TC89 


1988 Oct 14 


4.9 


< 1.08 


I&91 


1989 Oct 4 


8.3 


2.47 ±0.15 


S&93 


1994 Apr 7 


88.2 


< 7.1 


I&95 


1998 Feb 15 


4.9 


0.91 ±0.08 


Here 



ST90 - Seaquist & Taylor 1990; TC89 - Torbett & Campbell 
1989; I&91 - Ivison ct al. 1991; S&93 - Seaquist et al. 1993; I&95 
- Ivison et al. 1995. 

also changed: the intensity ratios of the [O ill] and [Ne ill] 
lines were consistent with a low-density photoionised region, 
n ~ a few x 10 6 cm~ 3 , T c ~ 10 000 - 20 000 k, while the 
high-ionization emission lines required a low-density, hot, 
n e g 10 7 cm~ 3 , T c ~ 10 5 K, mechanically- heated region. 
Kenyon et al. (1991), Mikolajewska & Kenyon (1992) associ- 
ated these changes with a transition from an optically thick 
to an optically thin accretion disc following the decline from 
the 1975 eruption and suggested that the high-ionization for- 
bidden line emission might be excited in an accretion disk 
corona that formed when the outer disk became optically 
thin. This transition occured within a year or so, and the 
[Fevn] and [Mgv] emission lines remained strong until at 
least 1994 (Mikolajewska & Holowacz 2000). Kenyon et al. 
(1991) estimated the spherical radius of this highly ionised 
region as R ~ 25 (d/2kpc) 2/3 AU, which is surprisingly close 
to the blackbody radius, R ~ 37 (T e /10 5 )" - 5 (d/2kpc)AU 
appropriate for the 4.9-GHz flux density reported by Tor- 
bett & Campbell (1989). If these connections between the 
radio emission and the spectral changes are not coincidental 
then the high-ionization [Fevn], [Nev] and Mg [v] lines and 
the radio emission might have a common origin. 

3.2 Binary-star models 

Another possible scenario involves an asymmetric, ionised 
region produced when a stellar wind is ionised by a hot com- 
panion, as in the STB model. The geometry of the ionised re- 
gion is governed by a single parameter, X (the ionization pa- 
rameter) , which depends on the red giant mass-loss rate, the 
binary separation, and the Lyman continuum luminosity of 
the ionizing source. In this model, knowledge of the turnover 
frequency, v t , can be used to estimate the binary separation, 
a, which can then be compared with independently-known 
values (Table 3). For a range in X covering two orders of 
magnitude, the binary separation (within a factor of 2) is 
(Seaquist & Taylor 1990): 




where S t is the optically thin flux near the turnover, and 
d is the distance. For v% = 27 GHz and St ~ 12.6 mjy, we 
find a ra( jio ~ 40(d/kpc) AU. A comparison with the binary 
separation derived from spectroscopic orbit solution, a sp ~ 
2.2 AU (Kenyon et al. 1991) shows that the observed turnover 
frequency over-estimates the binary separation by a factor 



© 0000 RAS, MNRAS 000, 000-000 



A test of popular models for symbiotic stars 5 



Table 3. Binary parameters of CI Cyg 





Adopted value 


Reference 


Distance [pc] 


2000 


K&91 


Orbital period [days] 


855.6 


B&00 


Mass ratio, M g /M h 


3 


K&91 


Inclination 


k, 73° 


K&91 


Separation of the stars [AU] 


2.2 


K&91 


Mass of the red giant [Mq] 


1.5 


K&91 


Mass of the hot component [Mq] 


0.5 


K&91 



K&91 - Kenyon et al. 1991; B&00 - Bclczyhski et al. 2000. 



of ~ 30 for d ~ 2 kpc. The STB model apparently does not 
account satisfactorily for the radio emission from CI Cyg. 

There are also other inconsistencies between the ob- 
served spectrum of CI Cyg and the STB model. We have 
already mentioned (see §3 above) the serious disagreement 
between the Lyman continuum luminosity derived from the 
optically thin radio emission and that estimated from the 
UV continuum and emission-line studies. In addition, the 
observed optically thick spectral index (a > +0.6) imme- 
diately implies X < -k/4 in the STB model. Since CI Cyg 
is a well-studied eclipsing system we know that our radio 
observations were made when the system was viewed along 
the binary axis with the hot star in front (i.e. the view- 
ing angle, as defined in STB, 6 « 0°). Thus a = +0.96 
indicates X w 0.7 (c.f. Fig. 3 of Seaquist & Taylor 1984). 
The nebula in CI Cyg should thus be radiation bounded. 
This is inconsistent with its observed forbidden-line spec- 
trum. In particular, after 1984 the only strong lines have 
been those with high excitation ([Fevn], [Nev] and [Mgv]) 
while the intermediate- and low-excitation lines have been 
very faint ([Oiii], [Nein]) or absent (e.g. [Oil], [Sn]) sug- 
gesting that the nebula is density bounded (X > 1). Further 
insight into the geometry of the nebular region(s) of CI Cyg 
is provided by the apparent lack of Raman-scattered O VI 
lines at AA6825, 7082. These lines, commonly observed in 
high-excitation symbiotic systems, have never been detected 
in CI Cyg, even at the epoch when the high-ionization lines 
were very strong. A plausible explanation seems to be a lack 
of sufficient neutral H° scatterers, which also implies a high 
value of X, 10). 

There are several plausible explanations for the failure 
of the STB model in the case of CI Cyg. Above all, CI Cyg 
is one of the few symbiotic binaries in which the cool giant 
fills or nearly fills its Roche lobe (e.g. Mikolajewska 1996, 
and references therein). This has two important implica- 
tions: first, the size of the giant fills a large fraction of the 
binary separation (~ 0.47a for q = M g /Mh ~ 3) and for 
X <, 1/3 the ionisation front impinges on its photosphere; 
second, and perhaps most important, mass loss from the gi- 
ant is strongly concentrated in a stream flowing through the 
inner langrangian point (Li) towards the hot component (as 
opposed to the spherically symmetric wind assumed in the 
STB model). The result is a complex density distribution, 
not the simple ~ r~ 2 distribution in the STB model. In 
particular, the density should vary roughly as r" 1 along the 
stream, possibly in addition to a more symmetric component 
with much lower density - an order of magnitude or more 
lower that that in the stream (e.g. Lubov & Shu 1975). Thus 



the simple geometry assumed for the STB model is incom- 
patible with the conditions in the circumbinary environment 
of CI Cyg. 

The fact that the giant component of CI Cyg fills its 
tidal lobe also implies that the interacting wind model (as- 
suming spherically symmetric winds from both components) 
cannot be a good match to the true conditions - the stream 
should result in the formation of a disk of material orbiting 
the companion. Any wind from the companion will thus be 
bipolar rather than spherically symmetric. 

Radio continuum spectra of collimated ionised stellar 
winds have been calculated by Reynolds (1986) showing that 
unresolved sources can have partially opaque spectra with 
spectral indices between +2 and —0.1. These models as- 
sume the thermal material appears at a finite radius, ro, 
where the jet half- width is too. Power-law dependencies of 
quantities on jet length, r, are also assumed; in particu- 
lar, the jet half-width, temperature, velocity, density and 
ionised fraction are assumed to vary as r/ro to powers of 
e, gx, (jV, g n , and g x , respectively. The models are defined 
by values of e, g v , and g x ; the other gradients are derived 
from these. In particular, a spectrum with a = +0.95, 
as observed for CI Cyg, can be produced by some of the 
models. A good example is a conical, recombining, cooling, 
constant-velocity jet: the observed spectral index can be re- 
produced by the following combination of gradients: g v = 0, 
e = 1, Ox = 0T = — 1/3 (using Reynolds' notation). Then 
q T — —3.22, and F(q T , a op ) = 1.24. Assuming the jet is per- 
pendicular to the orbital plane of CI Cyg (and the eclipses 
observed in CI Cyg do imply i ~ 90°) and using our ob- 
served values of the turnover frequency and the radio flux 
in equations (18) and (19) from Reynolds, we obtain a char- 
acteristic jet extension (collimation distance) and mass-loss 
rate of 

-~ 133 (nf'" (£)"'*»• < 5 > 

and 

M ~ 2 - 85Xl °"(2l)0^ < 6 > 

The jet would thus have an extension of 2ro/d ~ 133/42 
milliarcsec for T ~ 10 4 /10 5 K, respectively, and in principle 
should be detectable. 

It is interesting that for a reasonable O (~ 0.4, consis- 
tent with the ratio of the equatorial to spherical radius of 
H + and He +2 estimated from narrow eclipse features - see 
§3.1), the jet cuts out (0/2) 2 ~ 0.04 of the solid angle, com- 
parable to the ratio of the Lyman continuum luminosities 
as derived from the radio and UV data, respectively. 

According to Reynolds' models, the spectral index of 
the integrated flux reflects the gradients in the region of 
the flow corresponding to the frequency of observation; spa- 
tial structure in the source is thus reflected in the spectral 
signature of the flux. In general, the conditions may vary 
along the flow, resulting in steepening and flattening trends 
in spectra - an example is shown in Figure 2 of Reynolds 
(1986). If the inner jet of CI Cyg is confined (e < 1) it can 
be responsible for the flat spectrum (optically thin) between 
2 and 0.85 mm, whilst the cm range would be dominated by 
the free (e = 1) outer jet. 



© 0000 RAS, MNRAS 000, 000-000 



6 J. Mikolajewska & R. J. Ivison 



4 CONCLUSIONS 

We conclude that neither the STB nor the colliding winds 
model can account satisfactorily for the radio emission from 
CI Cyg. This is probably because the cool component in this 
binary fills its Roche lobe so that the geometry of the mass 
flow differs significantly from the simple geometries assumed 
for those models. 

It seems that the bulk of the Balmer H I emission is 
formed in a dense region between the binary components, 
probably in the partially ionised mass-loss stream, while 
the high-excitation forbidden lines and radio emission arises 
from lower density regions, perhaps in a bipolar outflow. 
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